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, Abstract. 

T-H ' In cosmography, cosmokinctics, and cosmology it is quite common to encounter 

■ physical quantities expanded as a Taylor series in the cosmological redshift z. Perhaps 
O ! the most well-known exemplar of this phenomenon is the Hubble relation between 
^ '' distance and redshift. However, we now have considerable high-z data available, for 
^ . instance we have supernova data at least back to redshift z « 1.75. This opens up 

I I ' the theoretical question as to whether or not the Hubble series (or more generally any 

series expansion based on the z-redshift) actually converges for large redshift? Based 

"kjJ ! on a combination of mathematical and physical reasoning, we argue that the radius of 

' convergence of any series expansion in z is less than or equal to 1, and that z-based 

OO . expansions must break down for z > 1, corresponding to a universe less than half its 
no . 

^ current size. 

• ' Furthermore, we shall argue on theoretical grounds for the utility of an improved 

, parameterization y — z/(l + z). In terms of the y-redshift we again argue that the 

' radius of convergence of any series expansion in y is less than or equal to 1, so that 

■ y-based expansions are likely to be good all the way back to the big bang {y = 1), but 
^ , that y-based expansions must break down for y < —I, now corresponding to a universe 

' more than twice its current size. 

X: 

^ , Keywords: high redshift, convergence. 
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1. Introduction 



Consider tlie standard luminosity distance versus redshift relation [1, 2]: 




(1) 



and its higher-order extension [3, 4, 5, 6] 



(2) 



As will quickly be verified by looking at the derivation (see, for example, [1, 2, 3, 4, 5, 6], 
the standard Hubble law is actually a Taylor series expansion derived for small z, whereas 
much of the most interesting recent supernova data occurs at 2 > 1 [7, 8, 9, 10, 11, 12]. 
Should we even trust the usual formalism for large z > 17 Two distinct things could go 
wrong [13]: 

• The underlying Taylor series could fail to converge. 

• Finite truncations of the Taylor series might be a bad approximation to the exact 



In fact, both things happen. There are good mathematical and physical reasons for this 
undesirable behaviour, as we shall discuss below. We shall carefully explain just what 
goes wrong — and suggest various ways of improving the situation. Our ultimate goal 
will be to find suitable forms of the Hubble relation that are well adapted to performing 
fits to all the available distance versus redshift data. More generally, the same sort of 
argument applies to any physical quantity that is expanded as a Taylor series in the 
redshift. Based on a combination of mathematical and physical reasoning, we argue 
that the radius of convergence of any series expansion in z is less than or equal to 1, 
and that 2-based expansions must break down for z > 1, corresponding to a universe 
less than half its current size. Furthermore, z-based expansions are likely to be good all 
the way to z = —1, corresponding to a universe that has grown to infinite size. 

We shall then argue on theoretical grounds for the utility of an improved 
parameterization y = z/{l + z). In terms of the redshift we again argue that the 
radius of convergence of any series expansion in y is less than or equal to 1, so that 
y-based expansions are likely to be good all the way back to the big bang {y = 1), but 
that |/-based expansions must break down for ?/ < — 1, now corresponding to a universe 
more than twice its current size. Choosing to adopt the 2;-redshift or the y-redshift will 
not change the physics, but it might improve mathematical convergence and physical 
insight. 

2. Cosmological distance scales 



result. 



In cosmology there are many different and equally natural definitions of the notion of 
"distance" between two objects or events, whether directly observable or not. For the 



The Hubble series: Convergence properties and redshift variables 



3 



vertical axis of the Hubble plot, instead of using the standard default choice of luminosity 
distance di, let us now consider using one or more of: 

• The "photon flux distance" : 

= {i + ly/^- 

• The "photon count distance" : 

• The "deceleration distance": 

<^Q = (i+%/2- (5) 

• The "angular diameter distance" : 

• The "distance modulus": 

fXD = 5 logioML/(10 pc)] = 5 logio[rfL/(l Mpc)] + 25. (7) 

• Or possibly some other surrogate for distance. 

Some words of explanation and caution are in order here [13, 14]: 

• The "photon flux distance" dp is based on the fact that it is often technologically 



easier to count the photon flux (photons/sec) than it is to bolometrically measure 
total energy flux (power) deposited in the detector. If we are counting photon 
number flux, rather than energy flux, then the photon number flux contains one 
fewer factor of (1 + z)~^. Converted to a distance estimator, the "photon flux 
distance" contains one extra factor of (1 + z)~^^'^ as compared to the (power-based) 
luminosity distance. 

The "photon count distance" dp is related to the total number of photons absorbed 
without regard to the rate at which they arrive. Thus the "photon count distance" 
contains one extra factor of (1 + z)~^ as compared to the (power-based) luminosity 
distance. Indeed D'Inverno [16] uses what is effectively this photon count distance as 
his nonstandard definition for luminosity distance. Furthermore, though motivated 
very differently, this quantity is equal to Weinberg's definition of proper motion 
distance [1], and is also equal to Peebles' version of angular diameter distance [2]. 
That is: 

dp '^L, D'Inverno '^proper, Weinberg '^A, Peebles- (8) 

The quantity dq is (as far as we can tell) a previously un-named quantity that seems 
to have no simple direct physical interpretation — but we shall soon see why it is 
potentially useful, and why it is useful to refer to it as the "deceleration distance" . 
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The quantity dA is Weinberg's definition of angular diameter distance [1], 
corresponding to tlie physical size of the object when the light was emitted, divided 
by its current angular diameter on the sky. This differs from Peebles' definition 
of angular diameter distance [2], which corresponds to what the size of the object 
would be at the current cosmological epoch if it had continued to co-move with 
the cosmological expansion (that is, the "comoving size"), divided by its current 
angular diameter on the sky. Weinberg's dA exhibits the (at first sight perplexing, 
but physically correct) feature that beyond a certain point dA can actually decrease 
as one moves to older objects that are clearly "further" away. In contrast Peebles' 
version of angular diameter distance is always increasing as one moves "further" 
away. Note that 

dA,Peehlcs = (1 + z) dA- (9) 

Finally, note that the distance modulus can be rewritten in terms of traditional 
stellar magnitudes as 

f^D /^apparent /^absolute- 

(10) 

The continued use of stellar magnitudes and the distance modulus in the context 
of cosmology is largely a matter of historical tradition, though we shall soon see 
that the logarithmic nature of the distance modulus has interesting and useful side 
effects. Note that we prefer as much as possible to deal with natural logarithms: 
Ina; = In(lO) log^g ^- Indeed 

fiD = T^ ln[rfi/(l Mpc)]+25, (11) 
m 10 



so that 



Obviously 



\n[dL/{l Mpc)] = - 25]. (12) 



dh ^ dp > dp > dq > dA- (13) 

Furthermore these particular distance scales satisfy the property that they converge on 
each other, and converge on the naive Euclidean notion of distance, as ^ — * 0. 

To simplify subsequent formulae, it is now useful to define the "Hubble distance" X 

du = ^, (14) 

SO that for Hq = 73 (km/sec)/Mpc [15] we have 

dn = 4100 tlfo Mpc- (15) 

I The "Hubble distance" dn = c/Hq is sometimes called the "Hubble radius", or the "Hubble sphere", 
or even the "speed of light sphere" [SLS] [17]. Sometimes "Hubble distance" is used to refer to the naive 
estimate d ~ dn z coming from the linear part of the Hubble relation and ignoring all higher-order 
terms — this is definitely not our intended meaning. 
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Furthermore we choose to set 

"« = l + W=l + ^- (16) 

For our purposes is a purely cosmographic definition without dynamical content. 
(Only if one additionally invokes the Einstein equations in the form of the Friedmann 
equations does Qq have the standard interpretation as the ratio of total density to the 
Hubble density, but we would be prejudging things by making such an identification in 
the current cosmographic framework.) In the cosmographic framework k/a^ is simply 
the present day curvature of space (not spacetime), while = Hq/c^ is a measure 
of the contribution of expansion to the spacetime curvature of the FLRW geometry. 
More precisely, in a FRLW universe the Riemann tensor has (up to symmetry) only two 
non-trivial components. In an orthonormal basis: 
_ k _ k 

Rm4> - ^ + ^- ^ + ^; (17) 

D '<^_qH^ ...^ 

Then at arbitrary times Q can be defined purely in terms of the Riemann tensor of the 
FLRW spacetime as 



3. New versions of the Hubble law 



New versions of the Hubble law are easily calculated for each of these cosmological 
distance scales. Explicitly: 



dd 


^) 


= dn z< 


1 - 




f ^ [go + 'iql - (jo + ^^o)] + 0(z=^)|. 


(20) 


dp 




= dn z< 


1 - 


1 1 r 
-^qoz + - [3 


+ lOgo + I2ql - 4(jo + l^o)] + 0(;2^)|. 


(21) 


dp 




= dn z< 


1 - 


^ [1 + 9o] ^ + 


^ [3 + 4go + H - (jo + fio)] z^ + 0{z^)^. 


(22) 


dq 


{^) 


= du Z' 


1 - 




^ [27 + 22go + I2ql - 4(jo + f^o)] z'' + 0{z^)^ 


. (23) 


dA 




= dn z< 


1 - 


^ [3 + go] ^ + 


i [12 + 7go + H - (jo + ^^o)] 2' + 0(^3)|, 


(24) 
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If one simply wants to deduce (for instance) the sign of Qq, then it seems that plotting 
the "photon flux distance" dp versus z would be a particularly good test — simply check 
if the first nonlinear term in the Hubble relation curves up or down. 

In contrast, the Hubble law for the distance modulus itself is given by the more 
complicated expression 

fioiz) = 25 + ^ I ln((ijy/Mpc) + In^; 



In(lO) 

+ l[l-qo]z-^[3- lOgo - 9ql + 4(jo + ^o)] + Oiz')^ (25) 

However, when plotting fi£, versus z, most of the observed curvature in the plot comes 
from the universal (In z) term, and so carries no real information and is relatively 
uninteresting. It is much better to rearrange the above as: 

ln[dL/{z Mpc)] = - 25] - In^ 

5 

= ln(c/jy/Mpc) 

- i [-1 + go] ^ + ^ [-3 + lOgo + 9ql - 4(jo + ^^o)] + 0{z'). (26) 
In a similar manner one has 

\n[dF/{z Mpc)] = i^[/iB-25]-ln^-^ln(l + z) 





\n{dH/Mpc 
- + 24 



l:qoZ + 7^ [3 + lOgo + 9ql - 4(jo + l^o)] + 0{z'). (27) 



ln[dp/{z Mpc)] = i^[/iD - 25] - In z - ln(l + z) 
5 

= ln(c/H/Mpc) 



^ [1 + go] 2 + ^ [9 + lOgo + 9go' - 4(jo + l^o)] 2' + 0{z^). (2^ 



\n[dQ/{z Mpc)] = i^[/i^-25]-ln2;-^ln(l + 2;) 
= ln(rfjy/Mpc) 



- ^ [2 + go] z + ^ [15 + lOgo + Dgo^ - 4(jo + ^o)] + 0{z^). (29) 
ln[ciA/(^ Mpc)] = i^[/iD - 25] - In ^ - 2 ln(l + z) 



= ln(dH/Mpc) 

- ^ [3 + go] 2 + ^ [21 + lOgo + 9g2 - 4(jo + ^o)] z^ + Oiz^). (30) 

These logarithmic versions of the Hubble law have several advantages — fits to these 
relations are easily calculated in terms of the observationally reported distance moduli 
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fin and their estimated statistical uncertainties [7, 8, 9, 10, 11]. (Specifically there is no 
need to transform the statistical uncertainties on the distance moduli beyond a universal 
multiplication by the factor [In 10]/5.) Furthermore the deceleration parameter go is easy 
to extract as it has been "untangled" from both Hubble parameter and the combination 

Note that it is always the combination (jo + ^o) that arises in these third-order 
Hubble relations, and that it is even in principle impossible to separately determine jo 
and Qq in a cosmographic framework. The reason for this degeneracy is (or should be) 
well-known [1, p. 451]: Consider the exact expression for the luminosity distance in any 
FLRW universe, which is usually presented in the form [1, 2] 

dL{z) = ao {1 + z) sink I I -^d^l, (31) 

Ho ao Jo H{z) J 

where 

{sin(x), k = +1; 
X, /c = 0; (32) 

sinh(a:), k = —1. 

By inspection, even if one knows H{z) exactly for all z one cannot determine diiz) 
without independent knowledge of k and ao. Conversely even if one knows diiz) exactly 
for all z one cannot determine H{z) without independent knowledge of k and oq. Indeed 
let us rewrite this exact result in a slightly different fashion as 

I ao Jo H{z) 

dL{z) = ao (1 + z) ^ -j= (33) 

where this result now holds for all k provided we interpret the /c = case in the obvious 
limiting fashion. Equivalently, using the cosmographic VLq as defined above we have the 
exact cosmographic result that for all Vlo'. 



di,(z) =d„(l + z) ^ ' ' (34) 

v"o ~ J- 

This form of the exact Hubble relation makes it clear that an independent determination 
of VLo (equivalently, k/a^), is needed to complete the link between a(t) and dL^z). When 
Taylor expanded in terms of z, this expression leads to a degeneracy at third-order, 
which is where f2o [equivalently k/a^ first enters into the Hubble series [5, 6]. 

What message should we take from this discussion? There are many physically 
equivalent versions of the Hubble law, corresponding to many slightly different physically 
reasonable definitions of distance, and whether we choose to present the Hubble law 
linearly or logarithmically. If one were to have arbitrarily small scatter /error bars on 
the observational data, then the choice of which Hubble law one chooses to fit to would 
not matter. In the presence of significant scatter/uncertainty there is a risk that the fit 
might depend strongly on the choice of Hubble law one chooses to work with. (And if 
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the resulting values of the deceleration parameter one obtains do depend significantly 
on which distance scale one uses, this is evidence that one should be very cautious 
in interpreting the results.) Note that the two versions of the Hubble law based on 
"photon flux distance" dp stand out in terms of making the deceleration parameter 
easy to visualize and extract. 

4. Why is the redshift expansion badly behaved for z > 11 

In addition to the question of which distance measure one chooses to use, there is a 
basic and fundamental physical and mathematical reason why the traditional redshift 
expansion breaks down for z > 1. 

4-1. Convergence 

Consider the exact Hubble relation (31). This is certainly nicely behaved, and possesses 
no obvious poles or singularities, (except possibly at a turnaround event where H{z) —>■ 
0, more on this below). However if we attempt to develop a Taylor series expansion in 
redshift z, using what amounts to the definition of the Hubble Hq, deceleration qq, and 
jerk jo parameters, then: 

T^ = ^ = l + Ho{t-to)-'-^{t- tor + ^ (t - to)' + 0{[t- to]% (35) 

Now this particular Taylor expansion manifestly has a pole at z = — 1, corresponding to 
the instant (either at finite or infinite time) when the universe has expanded to infinite 
volume, a = oo. Note that a negative value for z corresponds to a{t) > ao, that is: In 
an expanding universe z < corresponds to the future. Since there is an explicit pole 
at z = —1, by standard complex variable theory the radius of convergence is at most 
\z\ = 1, so that this series also fails to converge for z > 1, when the universe was less 
than half its current size. 

Consequently when reverting this power series to obtain lookback time T = to — ta.s 
a function T{z) of z, we should not expect that series to converge for 2; > 1. Ultimately, 
when written in terms of oq, Hq, qo, jo, and a power series expansion in redshift z you 
should not expect dii^z) to converge for 2 > 1. 

Note that the mathematics that goes into this result is that the radius of 
convergence of any power series is the distance to the closest singularity in the complex 
plane, while the relevant physics lies in the fact that on physical grounds we should not 
expect to be able to extrapolate forwards beyond a = 00, corresponding to z = — 1. 
Physically we should expect this argument to hold for any observable quantity when 
expressed as a function of redshift and Taylor expanded around z = — the radius of 
convergence of the Taylor series must be less than or equal to unity. (Note that the radius 
of convergence might actually be less than unity, this occurs if some other singularity 
in the complex z plane is closer than the breakdown in predictability associated with 
attempting to drive a(t) "past" infinite expansion, a = 00.) Figure 1 illustrates the 
radius of convergence in the complex plane of the Taylor series expansion in terms of z. 
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Complex z plane 


z = 


= -1 Z-- 


-- 


z - 


1= 1 z = ^ 


a 4 


F +00 a = 


= ao 


a i 




1= ao/2 a = 



radius of convergence 

Figure 1. Qualitative sketch of the behaviour of the scale factor a and the radius of 
convergence of the Taylor series in z-redshift. 

Consequently, we must conclude that observational data regarding dL^z) for 2; > 1 
is not going to be particularly useful in fitting oq, Hq, qq, and jo, to the usual traditional 
version of the Hubble relation. 



4-. 2. Pivoting 

A trick that is sometimes used to improve the behaviour of the Hubble law is to Taylor 
expand around some nonzero value of z, which might be called the "pivot" . That is, we 
take 

z = Zpi^ot + ^z, (36) 
and expand in powers of Az. If we choose to do so, then observe 

\. = 1 + ^0 it-to)-l go H'o it-tof + ^^ Jo Hi {t-tof + 0{[t-tof). (37) 
The pole is now located at: 

Az = -{1 + Zpi^ot). (38) 

which again physically corresponds to a universe that has undergone infinite expansion, 
a = 00. The radius of convergence is now 

\Az\<{l + Zpi,ot). (39) 

and we expect the pivoted version of the Hubble law to fail for 

z>l + 2 Zpiyot- (40) 

So pivoting is certainly helpful, and can in principle extend the convergent region of the 
Taylor expanded Hubble relation to somewhat higher values of z, but maybe we can do 
even better? 
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4-3. Other singularities 

Other singularities that might further restrict the radius of convergence of the Taylor 
expanded Hubble law (or any other Taylor expanded physical observable) are also 
important. Chief among them are the singularities (in the Taylor expansion) induced 
by turnaround events. If the universe has a minimum scale factor amin (corresponding 
to a "bounce") then clearly it is meaningless to expand beyond 

implying that we should restrict our attention to the region 

\z\ < 2;max = Cto/amm " 1- (42) 

Since for other reasons we had already decided we should restrict attention to \z\ < 1, 
and since on observational grounds we certainly expect any "bounce" , if it occurs at all, 
to occur for z^ax ^ 1, this condition provides no new information. 

On the other hand, if the universe has a moment of maximum expansion, and then 
begins to recoUapse, then it is meaningless to extrapolate beyond 

1 + Zmin = Oo/Omaxi Z^i^ = —[1 — do/o-max]; (43) 

implying that we should restrict our attention to the region 

< 1 - ao/Omax- (44) 

This relation now does provide us with additional constraint, though (compared to the 
1^1 < 1 condition) the bound is not appreciably tighter unless we are "close" to a point 
of maximum expansion. Other singularities could lead to additional constraints. 

5. Improved redshift variable for the Hubble relation 

Now it must be admitted that the traditional redshift has a particularly simple physical 
interpretation: 

l + . = ^ = 454, (45) 

Ae a[te) 

so that 

Ao — Ae AA 

z = — = — . (46) 

Ae 

That is, z is the change in wavelength divided by the emitted wavelength. This is 
certainly simple, but there's at least one other equally simple choice. Why not use: 
Ao - Ae AA 

y = — 7 = ^ ? (47) 

Ao Ao 

That is, define y to be the change in wavelength divided by the observed wavelength. 
This implies 

l-,= ^ = 4^4= 1 (48) 

Aq a{to) l + z 
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Now similar expansion variables have certainly been considered before. (See, for 
example, Chevalier and Polarski [18], who effectively worked with the dimensionless 
quantity b = a(t)/ao, so that y = 1 — b. Similar ideas have also appeared in several 
related works [19, 20, 21, 22]. Note that these authors have typically been interested 
in parameterizing the so-called w-parameter, rather than specifically addressing the 
Hubble relation.) 

Indeed, the variable y introduced above has some very nice properties: 

y = T^; z = T^- (49) 

1 + z 1 — 2/ 

In the past (of an expanding universe) 

^e(0,oo); ye (0,1); (50) 

while in the future 

ze(-l,0); ye (-00,0). (51) 

So the variable y is both easy to compute, and when extrapolating back to the Big Bang 
has a nice finite range (0, 1). We will refer to this variable as the y -redshift. (Originally 
when developing these ideas we had intended to use the variable y to develop orthogonal 
polynomial expansions on the finite interval y e [0,1]. This is certainly possible, but 
we shall soon see that given the current data, this is somewhat overkill, and simple 
polynomial fits in y are adequate for our purposes.) 

In terms of the variable y it is easy to extract a new version of the Hubble law by 
simple substitution: 

diiy) =dHy!^l-l[-S + qo]y + l [12 - 5qo + 3ql - (jo + ^^o)] + 0(y3)|. (52) 

This still looks rather messy, in fact as messy as before — one might justifiably ask in 
what sense is this new variable any real improvement? 

First, when expanded in terms of y, the formal radius of convergence covers much 
more of the physically interesting region. Consider: 

l-y=l+Hoit-to)-^qoH^ (t- t^f + ^ jo {t - tof + 0([t - to]')- (53) 

This expression now has no poles, so upon reversion of the series lookback time T = t^—t 
should be well behaved as a function T{y) oi y — at least all the way back to the Big 
Bang. (We now expect, on physical grounds, that the power series is likely to break 
down if one tries to extrapolate backwards through the Big Bang.) Based on this, we 
now expect d^i^y), as long as it is expressed as a Taylor series in the variable y, to be a 
well-behaved power series all the way to the Big Bang. In fact, since 

y = +1 ^ Big Bang, (54) 

we expect the radius of convergence to be given by \y\ = 1, so that the series converges 
for 

\y\ < 1- (55) 
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Consequently, when looking into the future, in terms of the variable y we expect to 
encounter problems at y = —1, when the universe has expanded to twice its current 
size. Figure 2 illustrates the radius of convergence in the complex plane of the Taylor 
series expansion in terms of y. 



y = -oo y =1 


= -1 y-- 


^ Compl 

= y-- 


ex y plane 
i= 1 


a = +00 a =S 


F 2ao a = ao a 


] ^ 

U 



radius of convergence 

Figure 2. Qualitative sketch of the behaviour of the scale factor a and the radius of 
convergence of the Taylor series in y-redshift. 

Note the tradeoff here — z is a useful expansion parameter for arbitrarily large 
universes, but breaks down for a universe half its current size or less; in contrast y is 
a useful expansion parameter all the way back to the Big Bang, but breaks down for 
a universe double its current size or more. Whether or not y is more suitable than z 
depends very much on what you are interested in doing. This is illustrated in Figures 1 
and 2. For the purposes of this article we are interested in high-redshift supernovae — 
and we want to probe rather early times — so it is definitely y that is more appropriate 
here. Indeed the furthest supernova for which we presently have both spectroscopic data 
and an estimate of the distance occurs at z = 1.755 [10], corresponding to y = 0.6370. 
Furthermore, using the variable y it is easier to plot very large redshift datapoints. 
For example, (though we shall not pursue this point in this article), the Cosmological 
Microwave Background is located at -Zcmb = 1088, which corresponds to ?/cmb = 0.999. 
This point is not "out of range" as it would be if one uses the variable z. 



6. More versions of the Hubble law 



In terms of this new redshift variable, the "linear in distance" Hubble relations are: 
ddy) = dHy^l-\[-^ + go] y + \ [12 - 5go + H - (jo + ^^o)] y^ + (56) 



dpiy) 



r 1 1 

dnyll-- [-2 + go] y + — [27 - 14go + 12go' - 4(jo + ^^c 



o)] y^ + 0(y^) .(57) 
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dp{y) = rfn yjl - ^ [-1 + go] 2/ + ^ [3 - 2go + 3go' - (jo + ^o)] + Oiy')^ 
dqiy) = 2/|l - |y + ^ [3 - 2go + Wo - 4(jo + f^o)] y' + 



(5^ 



(59) 



dA{y) = dn yjl - ^ [1 + go] y + ^ [go + 3go' - (jo + ^o)] y' + (60) 

Note that in terms of the y variable it is the "deceleration distance" dg that has the 
deceleration parameter go appearing in the simplest manner. Similarly, the "logarithmic 
in distance" Hubble relations are: 

\n[dL/{y Mpc)] = - 25] -\ny 

= ln(c?H/Mpc) 

- i [-3 + go] 2/ + ^ [21 - 2go + - 4(jo + l^o)] y' + Oiy'). (61) 

\n[dF/{y Mpc)] = i^[/i^-25]-lny + ^ln(l-y) 
= ln(ciH/Mpc) 

- ^ [-2 + go] 1/ + ^ [15 - 2go + 9go' - 4(jo + f^o)] + 0{y^). (62) 



\n[dp/{y Mpc)] = ^[/Wi? - 25] - In y + ln(l - y) 
= ln(rfiy/Mpc) 

- ^ [-1 + go] y + ^ [9 - 2go + 9go' - 4(jo + l^o)] y' + 0{y^). (63) 



HdQ/{y Mpc)] = ^[/iD - 25] - Iny + ^ ln(l - y) 
= ln(o?H/Mpc) 

- ^go 1/ + ^ [3 - 2go + 9g2 - 4(jo + l^o)] 1/' + 0{y^). (64) 



ln[rfA/(2/ Mpc)] = ^ [/iD - 25] - In y + 2 ln(l - y) 
= ln(c////Mpc) 

- ^ [1 + go] 2/ + ^ [-3 - 2go + 9g2 - 4(jo + l^o)] + 0{y'). (65) 
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Again note that the "logarithmic in distance" versions of the Hubble law are attractive in 
terms of maximizing the disentangling between Hubble distance, deceleration parameter, 
and jerk. Now having a selection of Hubble laws on hand, for various definitions of 
cosmological distance and redshift, we can start to confront the observational data to 
see what it is capable of telling us. This is the basis of the analysis in [13], which is well 
outside the scope of the current article. 

7. Conclusions 

Our main conclusions are threefold: 

• The use of the 2;-redshift for 2; > 1 is likely to lead to mathematical problems — 
specifically any Taylor series in z will be guaranteed to diverge for z > 1, and so 
finite truncations will be poor approximations to the underlying physical function. 
This is not all that early in the evolution of the universe — indeed many galaxies 
and supernovae are seen in the region z > 1, so one ignores this issue at one's peril. 

• The use of the ?/-redshift, where y = z/{l + z), is very much to be encouraged for 
z > 1 (corresponding to y > 1/2). Taylor series in the redshift are likely to be 
well behaved all the way back to the big bang (corresponding to y = 1). 

• By combining the notions of 2;-redshift, redshift, and the many reasonably 
standard notions of "cosmological distance" that have appeared in the literature, 
it is possible to extract many different versions of the Hubble law. Which version 
of the Hubble law one chooses to adopt for any specific purpose will depend on the 
specific question being addressed. For instance, in [13] we have argued that plotting 
\n[dF/z Mpc] versus z, or \n[dQ/y Mpc] versus y, leads to a good visual test of the 
quality and robustness of the data underlying determinations of the deceleration 
parameter go- 
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